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Abstract 

Vibronic coupling constants in the monoanionic, trianionic, and excited states of picene are 



1-^ , evaluated from the total energy gradients using the density functional theory. Employing the 

calculated vibronic coupling constants in the excited state of the neutral molecule, electron energy 
^ ■ loss spectrum (EELS) is simulated to be compared with the experimental spectrum. The calculated 

o 
o 



vibronic coupling constants are analyzed in terms of the vibronic coupling density which enables 
us to analyze vibronic couplings based on the relation between the electronic and vibrational 
structures. The vibronic coupling constants reported by Kato et al. [J. Chem. Phys. 116, 
3420(2002) and Phys. Rev. Lett. 107, 077001 (2011)] are critically discussed based on the vibronic 
coupling density analysis. 

PACS numbers: 71.38.-k,79.20.Uv,33.20.Wr 



After the discovery of the superconductivity in alkah-metal (K, Rb) doped picene,- ex- 
perimental and theoretical studies on the electronic structure on the picene have been piled 
up.-""— The vibronic coupling (electron-vibration coupling)^ is an important interaction 
in the electronic properties such as superconductivity. Okazaki et al. have discussed an 
importance of the vibronic couplings in doped picene based on their photoelectron spec- 
tra (PES).- Therefore, evaluation of the vibronic coupling constants (VCC) is crucial to 
discuss electronic properties of doped picenes. Some authors have published the VCCs of 
picene anions.— >ii>i^'i^ However, some of the calculated VCCs are controversial. Subedi and 
Boeri have concluded that the electron-phonon coupling of the modes around 1600 cm~^ are 
strong,-lfi while those calculated by Kato et al are weak in this region.— li^ii^ 

Vibronic effects can be experimentally observed in spectra.— Roth et al. measured elec- 
tron energy loss spectrum (EELS) of pristine picene at 20 K.- They have observed vibronic 
progressions in the EELS of the intramolecular excitation to 5*2 (^-62) state. 

We have recently published calculation of the VCCs in CgQ from the gradients of the total 
energies.— The results are consistent with the experimental observation of the PES by Wang 
et alr^ We have proposed a concept, vibronic coupling density (VCD).— 1^ Based on the VCD, 
we can discuss vibronic couplings from view of electronic and vibrational structures. On the 
basis of the VCD analysis, we have succeeded in designing carrier-transporting molecules 
with small vibronic couplings which is required in organic electronics such as organic light- 
emitting diodes (OLED).-^'23 

In this work, we report the VCCs of the excited state of the free molecule in the neutral 
state, the monoanionic, and the trianionic states of the free molecule based on the same 
method of calculation employed in the calculation for C^q. Using the VCCs of the excited 
^i?2 state, we simulate EELS and compare the spectrum with the experimental oneP- 
the view of the VCD analysis, we critically discuss the previous VCCs in Refs. 
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We evaluated VCCs of mode a Va from the gradients of the adiabatic potential energy 
surface E with respect to a mass- weighted normal coordinates Qa'r^'^— 
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where H denotes a molecular Hamiltonian, Rq is the equilibrium geometry of the ground 
state of the neutral picene, \E' is an electronic wavefunction of the excited or anionic state at 



Ro, and the phase of a normal mode a is chosen so that Va becomes negative. The vibronic 
Hamiltonian is written as 

1 



-nvibro / J 



T{Qa) + 2^'Q' + Vo^Qc 



(2) 



where T{Qa) denotes kinetic energy operator of a vibration and Ua vibrational frequency. 
We employed Becke's hybrid functional (B3LYP)^ and Perdew and Wang's one with gen- 
eralized gradient approximation (PW91)^ with the triple-zeta 6-311+G(d,p) basis set. The 
geometries were optimized for the neutral ground state. The optimized structures with C2v 
symmetry were checked with vibrational analysis to be a minimum. The time-dependent 
density-functional-theory is applied for the excited state calculations. We performed ana- 
lytical force-calculations for the excited and the anionic states to obtain the VCCs.— The 
electronic and vibrational structures as well as the forces were obtained using a program 
package, Gaussian 09.— The VCCs are calculated using our codes. 

Since the fine structure observed in the EELS by Roth et at— can originate from the 
vibronic couplings which they did not take into account in Ref. |3|, we simulated the EELS 
to find an appropriate functional for the VCC calculations. EELS was simulated employing 



the same method as described in Ref. 118l We considered thermal excitation at 20K where 
Roth et al. observed the EELS. The calculated excitation energies are 3.7040 and 3.2520 eV 
for the B3LYP and PW91 functionals, respectively. The result using the PW91 functional 
reproduces well the experimental excitation energy 3.25 eV. Using the calculated VCCs in 
the 5*2(^-62) state, we simulated the EELS (FIGs. [1]). In the simulations, the 0-0 transition 
is set to 3.24 eV, and the linewidth a is assumed to be 270 cm~^ (The FWHM is 39.4 meV). 
The calculated VCCs^ as well as the vibrational frequencies are tabulated in Supporting 
Information. The spectrum calculated employing the PW91 functional (FIG.[l]^a)) shows a 
better fit than that using the B3LYP functional (FIG. [T](b)). In the calculation employing 
the B3LYP, the second strongest vibronic couplings around 1350 cm~^ are estimated larger 
than those in the result using the PW91.— Therefore, the relative intensities of EELS in the 
lower energy are reproduced employing the PW91.— Judging from the calculations of the 
excitation energies and the simulated spectra, we employ the PW91 functional hereafter for 
the calculation of the ionic states. The further comparisons of functionals are out of the 
scope in the present work. 

We calculated the VCCs V^ in the monoanion (^^2), dianion (^v4i), and trianion (^-Bi) 
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FIG. 1. (Color online) Red lines are simulated electron energy loss spectra (EELS) calculated by 
(a) the PW91 and (b) the B3LYP functional. Blue dots indicate the experimental one by Roth 
et a/.— The 0-0 transition is set to 3.24 eV. The other parameters employed in the simulations are 
T = 20 K and cr = 270 cm^^ 

(see Supporting Information). Kato et al. have calculated orbital vibronic coupling con- 
stants (OVCC) of the lowest unoccupied molecular orbitals (LUMO) and next LUMO 
(NLUMO).^^^ Hence we also calculated the OVCCs Vi^^ from the VCCs K. for com- 
parison with the ones previously reported. The vibronic Hamiltonian ([2]) is mapped onto 



a model Hamiltonian. The model Hamiltonian considered in Refs. 
follows: 
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is written as 



H. 
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(3) 



where orbitals i are LUMO (L) and NLUMO (NL), and a runs over all the active ai modes.— 
The dimensionless OVCC gi^a is defined by gi^a = ^,q/\/^^- ^L(^io-) is the creation (anni- 
hilation) operator of orbital i and spin o", and h]^{ha) is the creation (annihilation) operator 
of mode a. The OVCCs of LUMO 14,^ and NLUMO ^NL.a are obtained from the VCCs 
of monoanion Knono.a and the difference between VCCs of trianion Vtri,a and dianion Vdi^a, 
respectively: Vl,^ = Knono,a and Vnl.q = Vtn,a — Vdi,a. It should be noted that the present 
OVCCs effectively incorporate the contributions from all the occupied orbitals which are im- 
portant in quantitative arguments.— The OVCCs Vi^a and intramolecular electron-phonon 
couplings Aja/A^(0) = V^^/ul^ are shown in FIGs. [2] and tabulated in Supporting Informa- 
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FIG. 2. (Color online) Calculated (top) orbital vibronic coupling constants Vi^a and (bottom) 
electron-phonon couplings Aj^q/A^(0) in the (a) monoanion and the (b) trianion. The subscripts 
L and NL denote LUMO and NLUMO, respectively. Insets show the orbital vibronic couplings 
between 1300 to 1650 cni~^. 



tion. A^(0) is the density of states at the Fermi leveL 

In both the monoanion and trianion, two sets of modes with strong VCCs are observed: 
(1) the maximal coupling modes around 1600 cm~^ and (2) the second group of strong 
modes around 1350 cm~^. This is qualitatively consistent with the calculations by Subedi 
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IGL and 



17,— the modes around 1600 



and Boeri.— "^ However, in the calculations of Refs. 
cm~^ are weak for both the monoanion and trianion. The total electron-phonon couplings 
Al/A^(0) = Ea^L,a/iV(0) and ANL/iV(0) = Ea ^NL,a/A^(0) are 196.6 meV and 119.9 meV, 
respectively. The total coupling of the trianion Anl/^(0) is in line with that of Subedi 



and Boeri (110 ± 5 meV).^^ On the other hand, in Ref. [ll|, Al/A^(0) = 178 meV and 
-^nl/^IO) = 206 meV. They underestimated and overestimated the total couphngs for the 
monoanion and trianion, respectively. 

We will discuss the disagreement from view of the electronic and vibrational structures. 
The calculated VCCs can be rationalized based on the vibronic coupling density.—"^ A VCD 
r]a is defined as 

?7a(r) = Ap(r) x Va{r), (4) 

where Ap(r) = p(r) — po(r) is the electron density difference between the electron density of 
an ionic state p and that of a neutral state po. The potential derivative Va{r) is the derivative 
with respect to a mass- weighted normal coordinate Qa of the potential u{r) acting on one 
electron at a position r from all the nuclei. The vibronic coupling constant is equal to the 
integral of ?7a(r) over space r: 

K= [dhr]^{r). (5) 



The VCD gives a local picture of the vibronic coupling, and hence enables us to discuss the 
strength of the coupling qualitatively. 

We will concentrate on the ai(27) mode of the maximal-coupling mode around 1600 cm~^ 
and the ai(21) mode from the second group around 1350 cm~^ because the a-\ (27) and ai(21) 
modes are close to the mode 5 in Ref. [lO| and the 21st mode {1^21) in Ref. [ll|, respectively. 



FIGs. [3] show the potential derivatives Va{r) for the ai(27) and the ai(21) modes. The 
distribution of the potential derivative with respect to the ai(27) mode is located on the 
armchair-edges of the central three hexagons. On the other hand, that with respect to the 
ai(21) mode is on the terminal hexagons. 

FIG. 111(b) shows the electron density difference Ap of the monoanion. Since an additional 
electron occupies the LUMO (FIG. IH^a)), the positive vr density (white) appears in Ap. It 
should be noted that there occurs decrease of the a density (blue) in the molecular plane. 
Such a polarized density originates from the Coulomb interactions between the electron oc- 
cupying the LUMO and all electrons in doubly occupied orbitals below the highest occupied 
molecular orbital (HOMO).-^ 

FIGs. Oshow the vibronic coupling densities ^^^(r) of the monoanion with respect to the 
ai(27) and ai(21) modes. As for the ai(27), the electron density difference Ap(r) shows 
considerable overlap with the potential derivative Va{r). On the other hand, Ap(r) does not 
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FIG. 3. (Color online) Potential derivatives Va for (a) ai(27) (1605 cm ^) and (b) ai(21) (1379 
cni^^) modes. White and blue areas indicate positive and negative, respectively. The threshold is 
1.0 X 10-2 a.u. 
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FIG. 4. (Color online) (a) LUMO and (b)electron density difference Ap for the monoanion. White 
and blue areas indicate positive and negative, respectively. The threshold is 5.0 x 10^^ a.u. for 
LUMO and 1.8 x 10"^ a.u. for Ap. 

significantly overlap Wa(r) for the ai(21) mode. Therefore, the VCD rja (Hj) for the ai(27) is 
larger than that for the ai(21). Particularly, rja for the ai(27) has the large distribution on 
the central three armchair-edges. Accordingly, the VCC of the ai(27) mode is larger than 
that of the ai(21) mode. 

The VCD rja for the ai(27) mode on the bonds appears due to the polarization of Ap, thus 
neglecting such a polarization can give rise to quantitative, or sometimes qualitative, errors 
in VCC calculations based on the orbital levels.— In addition, an electron density difference is 
usually different from the orbital density of HOMO or LUMO. Many-body effect sometimes 



(a) ai(27) (1605 cm-i) 



(b) ai(21) (1379 cm-i) 




FIG. 5. (Color online) Vibronic coupling densities for (a) ai(27) (1605 cm ^) and (b) ai(21) (1379 
cm~^) modes of the monoanion. White and blue areas indicate positive and negative, respectively. 
The threshold is 2.0 x lO'^ a.u. 

TABLE I. Calculated frequencies (cm~^) and orbital gradients de/dQ and total-energy gradients 
dE/dQ as vibronic coupling constants (10~^ a.u.) for the selected modes of the monoanion. The 
subscripts L and S denote the neutral LUMO and the SOMO of the anion, respectively. 



Preq. dei^/dQ des/dQ dE/dQ 

Num. Anal. 

ai(2) 254.5 0.372 0.475 0.426 0.426 

ai(19) 1336.4 2.019 1.889 1.988 1.990 

ai(20) 1366.1 1.872 1.365 1.639 1.624 

ai(21) 1378.8 0.820 1.335 1.130 1.123 

ai(27) 1604.9 3.070 2.281 2.706 2.682 

ai(28) 1607.9 2.090 2.580 2.377 2.374 



plays a crucial role on the electron density difference, and therefore vibronic couplings.— 

Kato et al. have calculated VCCs as the gradients of the orbital levels (orbital gradients) 
with respect to normal coordinates.— We also obtained the orbital gradients of selected 
modes from the neutral LUMO dei^/dQa and from the singly occupied molecular orbital 
(SOMO) of the monoanion des/dQa to compare the present VCCs calculated analytically 
from the total energy gradient as well as numerical gradients. The orbital gradients and 

8 



the numerical energy gradients were obtained by fitting linear and quadratic polynomials, 
respectively, in the range from —0.2 to 0.2 a.m.u.^^^ ao with a step size 0.05 a.m.u.^^^ uq 
where Oq is the Bohr radius. These range and step size could be different from those of Kato 
et al. We summarize the results of the calculations in TABLE [H All the results indicate that 
the vibrational modes around 1600 cm~^ have the maximal coupling. 

Though the gradients of the LUMO level can yield results which are qualitatively consis- 
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and 



16 are not the case. 



tent with the gradients of the total energy, the results in Refs. 
In their calculations, the vibrational mode around 1380 cm~^ has the maximal coupling, 
2.3x10"^ a.u. for the monoanion.— However, since Ap (see FIGJl](b)) is mainly located 
on the armchair-edges of the central three hexagons, it does not overlap with the potential 
derivatives Va of the mode (FIG|3](b)), the VCC of the mode cannot be the maximal. Similar 
discussion holds for the mode around 1520 cm~^ which has the maximal coupling, 2.6x10"^ 
a.u. for the trianion.— 

In summary, we calculated the VCCs of picene for the excited state ^B2 of the neutral 
molecule, the monoanionic, and the trianionic states of the molecule from the gradients of 
the total energies with respect to the normal modes. In the previous studies, the VCCs 
calculated form the gradients of the orbital energies of the frontier level, LUMO. In other 
words, they have regarded the OVCCs of the LUMO as the VCCs. The present approach 
can provide quantitatively reliable VCCs, since all the occupied orbitals can contribute to 
the VCCs due to the selection rule. This has been discussed in Ref. Il8| in detail. We 
simulated the EELS and compared the spectra with experimental spectrum by Roth et al. 
No comparison with experiments has been reported on the vibronic couplings in picene. 
Needless to say, such a comparison between the theoretical and the experimental results is 
necessary to obtain reliable VCCs. From these simulations, we determined an appropriate 
functional for the calculation of the vibronic couplings in picene. Using the functional which 
reproduces the EELS, we evaluated the vibronic coupling constants of the picene anions. 
The calculated vibronic couplings of the anions can be employed for assignments of spectra, 
for example, photoelectron spectra of the anions. We discussed the vibronic couplings of the 
picene anions in terms of the vibronic coupling density. Based on the analysis, the present 
vibronic couplings are reasonable compared with the values reported previously. 
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